The genetic legacy of the expansion of Bantu-speaking peoples in Africa.

The expansion of people speaking Bantu languages is the most dramatic demographic event in Late Holocene Africa and fundamentally reshaped the linguistic, cultural and biological landscape of the continent1-7. With a comprehensive genomic dataset, including newly generated data of modern-day and ancient DNA from previously unsampled regions in Africa, we contribute insights into this expansion that started 6,000-4,000 years ago in western Africa. We genotyped 1,763 participants, including 1,526 Bantu speakers from 147 populations across 14 African countries, and generated whole-genome sequences from 12 Late Iron Age individuals8. We show that genetic diversity amongst Bantu-speaking populations declines with distance from western Africa, with current-day Zambia and the Democratic Republic of Congo as possible crossroads of interaction. Using spatially explicit methods9 and correlating genetic, linguistic and geographical data, we provide cross-disciplinary support for a serial-founder migration model. We further show that Bantu speakers received significant gene flow from local groups in regions they expanded into. Our genetic dataset provides an exhaustive modern-day African comparative dataset for ancient DNA studies10 and will be important to a wide range of disciplines from science and humanities, as well as to the medical sector studying human genetic variation and health in African and African-descendant populations.

Semipermeable species boundaries create opportunities for gene flow and adaptive potential.

Hybridisation and gene flow can have both deleterious and adaptive consequences for natural populations and species. To better understand the extent of hybridisation in nature and the balance between its beneficial and deleterious outcomes in a changing environment, information on naturally hybridising nonmodel organisms is needed. This requires the characterisation of the structure and extent of natural hybrid zones. Here, we study natural populations of five keystone mound-building wood ant species in the Formica rufa group across Finland. No genomic studies across the species group exist, and the extent of hybridisation and genomic differentiation in sympatry is unknown. Combining genome-wide and morphological data, we demonstrate more extensive hybridisation than was previously detected between all five species in Finland. Specifically, we reveal a mosaic hybrid zone between Formica aquilonia, F. rufa and F. polyctena, comprising further generation hybrid populations. Despite this, we find that F. rufa, F. aquilonia, F. lugubris and F. pratensis form distinct gene pools in Finland. We also find that hybrids occupy warmer microhabitats than the nonadmixed populations of cold-adapted F. aquilonia, and suggest that warm winters and springs, in particular, may benefit hybrids over F. aquilonia, the most abundant F. rufa group species in Finland. In summary, our results indicate that extensive hybridisation may create adaptive potential that could promote wood ant persistence in a changing climate. Additionally, they highlight the potentially significant ecological and evolutionary consequences of extensive mosaic hybrid zones, within which independent hybrid populations face an array of ecological and intrinsic selection pressures.

Genomic assessment reveals signal of adaptive selection in populations of the Spotted rose snapper Lutjanus guttatus from the Tropical Eastern Pacific.

Background: The lack of barriers in the marine environment has promoted the idea of panmixia in marine organisms. However, oceanographic conditions and habitat characteristics have recently been linked to genetic structure in marine species. The Tropical Eastern Pacific (TEP) is characterized by dynamic current systems and heterogeneous oceanographic conditions. The Gulf of Panama (part of the equatorial segment for the TEP) is influenced by a complex current system and heterogeneous environment, which has been shown to limit the gene flow for shoreline species. Next Generation Sequencing (NGS) has contributed to detect genetic differences in previously reported panmictic species by the assessment of loci associated with selection and to understand how selection acts affects marine populations. Lutjanus guttatus is a species distributed in the TEP for which previous studies using mitochondrial data recovered a panmictic pattern along its distributional range. In this study, we used SNP data of L. guttatus individuals sampled along its range to evaluate population genetic structure and investigate whether oceanographic factors influence the species' genetic architecture. Finally, we assessed the role of adaptive selection by evaluating the contribution of outlier and neutral loci to genetic divergence. Methods: The RADcap method was used to obtain 24 million paired reads for 123 individuals of L. guttatus covering nearly all its distributional area. Genetic variation was assessed using both spatial and non-spatial methods by comparing three different data sets: (i) a Combined Loci (CL dataset = 2003 SNPs); a search for putative loci under selection allowed the evaluation of (ii) Neutral Loci (NL dataset = 1858 SNPs) and (iii) Outlier Loci (OL dataset = 145 SNPs). We used the estimating effective migration surface (EEMS) approach to detect possible barriers to gene flow. Results: Genetic differences were found in the OL dataset, showing two clusters (Northern and Southern), whereas NL showed no differences. This result may be related to the Selection-Migration balance model. The limit between the Northern and Southern groups was in the Gulf of Panama, which has been previously identified as a barrier to gene flow for other species, mainly due to its heterogeneous oceanographic conditions. The results suggest that selection plays an important role in generating genetic differences in Lutjanus guttatus. A migration corridor was detected that coincides with the Costa Rica Coastal Current that flows from Central America to the Gulf of California, allowing the homogenization of the northern population. In the Southern cluster, a migration corridor was observed with the OL from Panama to Colombia, which could be associated with the currents found in the Gulf of Panama. Genetic variation found in the OL of Lutjanus guttatus highlights the usefulness of NGS data in evaluating the role of selection in population differentiation.

Host-plant adaptation as a driver of incipient speciation in the fall armyworm (Spodoptera frugiperda).

BACKGROUND: Divergent selection on host-plants is one of the main evolutionary forces driving ecological speciation in phytophagous insects. The ecological speciation might be challenging in the presence of gene flow and assortative mating because the direction of divergence is not necessarily the same between ecological selection (through host-plant adaptation) and assortative mating. The fall armyworm (FAW), a major lepidopteran pest species, is composed of two sympatric strains, corn and rice strains, named after two of their preferred host-plants. These two strains have been hypothesized to undergo incipient speciation, based on (i) several lines of evidence encompassing both pre- and post-zygotic reproductive isolation, and (ii) the presence of a substantial level of genetic differentiation. Even though the status of these two strains has been established a long time ago, it is still yet to be found whether these two strains indeed exhibit a marked level of genetic differentiation from a large number of genomic loci. Here, we analyzed whole genome sequences from 56 FAW individuals either collected from pasture grasses (a part of the favored host range of the rice strain) or corn to assess the role of host-plant adaptation in incipient speciation. RESULTS: Principal component analysis of whole genome data shows that the pattern of divergence in the fall armyworm is predominantly explained by the genetic differentiation associated with host-plants. The level of genetic differentiation between corn and rice strains is particularly marked in the Z chromosome. We identified one autosomal locus and two Z chromosome loci targeted by selective sweeps specific to rice strain and corn strain, respectively. The autosomal locus has both increased DXY and FST while the Z chromosome loci had decreased DXY and increased FST. CONCLUSION: These results show that the FAW population structure is dominated by the genetic differentiation between corn and rice strains. This differentiation involves divergent selection targeting at least three loci, which include a locus potentially causing reproductive isolation. Taken together, these results suggest the evolutionary scenario that host-plant speciation is a driver of incipient speciation in the fall armyworm.

Complex population history affects admixture analyses in nine-spined sticklebacks.

Introgressive hybridization is an important process in evolution but challenging to identify, undermining the efforts to understand its role and significance. On the contrary, many analytical methods assume direct descent from a single common ancestor, and admixture among populations can violate their assumptions and lead to seriously biased results. A detailed analysis of 888 whole-genome sequences of nine-spined sticklebacks (Pungitius pungitius) revealed a complex pattern of population ancestry involving multiple waves of gene flow and introgression across northern Europe. The two recognized lineages were found to have drastically different histories, and their secondary contact zone was wider than anticipated, displaying a smooth gradient of foreign ancestry with some curious deviations from the expected pattern. Interestingly, the freshwater isolates provided peeks into the past and helped to understand the intermediate states of evolutionary processes. Our analyses and findings paint a detailed picture of the complex colonization history of northern Europe and provide backdrop against which introgression and its role in evolution can be investigated. However, they also expose the challenges in analyses of admixed populations and demonstrate how hidden admixture and colonization history misleads the estimation of admixture proportions and population split times.

Evolutionary and ecological drivers of local adaptation and speciation in a North American avian species complex.

Throughout the speciation process, genomic divergence can be differentially impacted by selective pressures, as well as gene flow and genetic drift. Disentangling the effects of these evolutionary mechanisms remains challenging, especially for nonmodel organisms. Accounting for complex evolutionary histories and contemporary population structure often requires sufficient sample sizes, for which the expense of full genomes remains prohibitive. Here, we demonstrate the utility of partial-genome sequence data for range-wide samples to shed light into the divergence process of two closely related ducks, the Mexican duck (Anas diazi) and mallard (A. platyrhynchos). We determine the role of selective and neutral processes during speciation of Mexican ducks by integrating evolutionary and demographic modelling with genotype-environment and genotype-phenotype association testing. First, evolutionary models and demographic analyses support the hypothesis that Mexican ducks originally diverged ~300,000 years ago in climate refugia arising during a glacial period in southwest North America, and that subsequent environmental selective pressures played a key role in divergence. Mexican ducks then showed cyclical demographic patterns that probably reflected repeated range expansions and contractions, along with bouts of gene flow with mallards during glacial cycles. Finally, we provide evidence that sexual selection acted on several phenotypic traits as a co-evolutionary process, facilitating the development of reproductive barriers that initially arose due to strong ecological selection. More broadly, this work reveals that the genomic and phenotypic patterns observed across species complexes are the result of myriad factors that contribute in dynamic ways to the evolutionary trajectories of a lineage.

Ancient DNA and deep population structure in sub-Saharan African foragers.

Multiple lines of genetic and archaeological evidence suggest that there were major demographic changes in the terminal Late Pleistocene epoch and early Holocene epoch of sub-Saharan Africa1-4. Inferences about this period are challenging to make because demographic shifts in the past 5,000 years have obscured the structures of more ancient populations3,5. Here we present genome-wide ancient DNA data for six individuals from eastern and south-central Africa spanning the past approximately 18,000 years (doubling the time depth of sub-Saharan African ancient DNA), increase the data quality for 15 previously published ancient individuals and analyse these alongside data from 13 other published ancient individuals. The ancestry of the individuals in our study area can be modelled as a geographically structured mixture of three highly divergent source populations, probably reflecting Pleistocene interactions around 80-20 thousand years ago, including deeply diverged eastern and southern African lineages, plus a previously unappreciated ubiquitous distribution of ancestry that occurs in highest proportion today in central African rainforest hunter-gatherers. Once established, this structure remained highly stable, with limited long-range gene flow. These results provide a new line of genetic evidence in support of hypotheses that have emerged from archaeological analyses but remain contested, suggesting increasing regionalization at the end of the Pleistocene epoch.

Opportunities for studying propagule pressure using gene flow reveal its role in accelerating biological invasions.

When an alien species establishes at a new location, it must spread to become an invader. The extent to which propagule pressure promotes the spread of invaders, especially at local scales, is often difficult to quantify because it requires a reliable measure of, and variation in, rate of spread, and of propagule pressure across similar areas. In this issue of Molecular Ecology, Mairal et al. (2022) make use an unique system of paired sub-Antarctic islands, one with very infrequent human activities, and another inhabited by scientists, to assess the role of propagule pressure and anthropogenic disturbance in the introduction and spread of a major global invader, Poa annua L., to and on the islands. Genetic admixture between different genetic clusters is virtually absent from the little-visited island, while the inhabited island experienced more introduction events, but also significant admixture between genetic clusters. Detailed distribution maps of P. annua spanning more than 50 years allowed the authors to link genetic diversity to residence time. The nature of the system, and the multifaceted approach used by the authors, allows for new insights into the mechanism by which propagule pressure results in the spread of invasive species.

Ghost Lineages Highly Influence the Interpretation of Introgression Tests.

Most species are extinct, those that are not are often unknown. Sequenced and sampled species are often a minority of known ones. Past evolutionary events involving horizontal gene flow, such as horizontal gene transfer, hybridization, introgression, and admixture, are therefore likely to involve "ghosts," that is extinct, unknown, or unsampled lineages. The existence of these ghost lineages is widely acknowledged, but their possible impact on the detection of gene flow and on the identification of the species involved is largely overlooked. It is generally considered as a possible source of error that, with reasonable approximation, can be ignored. We explore the possible influence of absent species on an evolutionary study by quantifying the effect of ghost lineages on introgression as detected by the popular D-statistic method. We show from simulated data that under certain frequently encountered conditions, the donors and recipients of horizontal gene flow can be wrongly identified if ghost lineages are not taken into account. In particular, having a distant outgroup, which is usually recommended, leads to an increase in the error probability and to false interpretations in most cases. We conclude that introgression from ghost lineages should be systematically considered as an alternative possible, even probable, scenario. [ABBA-BABA; D-statistic; gene flow; ghost lineage; introgression; simulation.].

Repeated origins, widespread gene flow, and allelic interactions of target-site herbicide resistance mutations.

Causal mutations and their frequency in agricultural fields are well-characterized for herbicide resistance. However, we still lack understanding of their evolutionary history: the extent of parallelism in the origins of target-site resistance (TSR), how long these mutations persist, how quickly they spread, and allelic interactions that mediate their selective advantage. We addressed these questions with genomic data from 19 agricultural populations of common waterhemp (Amaranthus tuberculatus), which we show to have undergone a massive expansion over the past century, with a contemporary effective population size estimate of 8 x 107. We found variation at seven characterized TSR loci, two of which had multiple amino acid substitutions, and three of which were common. These three common resistance variants show extreme parallelism in their mutational origins, with gene flow having shaped their distribution across the landscape. Allele age estimates supported a strong role of adaptation from de novo mutations, with a median age of 30 suggesting that most resistance alleles arose soon after the onset of herbicide use. However, resistant lineages varied in both their age and evidence for selection over two different timescales, implying considerable heterogeneity in the forces that govern their persistence. Two such forces are intra- and inter-locus allelic interactions; we report a signal of extended haplotype competition between two common TSR alleles, and extreme linkage with genome-wide alleles with known functions in resistance adaptation. Together, this work reveals a remarkable example of spatial parallel evolution in a metapopulation, with important implications for the management of herbicide resistance.

Do we have infraspecific taxa of Salvia macrosiphon Boiss. (Lamiaceae) in Iran?

BACKGROUND: Salvia macrosiphon is an aromatic perennial species of Lamiaceae family that grows naturally in different parts of Iran. This herb is widely used in folk and modern medicine. Although in Flora Iranica and Flora of Iran, no infraspecific taxonomic rank has been detected for S. macrosiphon, some infraspecific taxonomic ranks have been reported. In the current study, we evaluated the genetic diversity and structure of 11 populations of this species to detect inter and intrapopulation genetic diversity and to survey the possibility of infraspecific taxonomic ranks in this species. METHODS AND RESULTS: We utilized the modified C-TAB protocol for DNA extraction and amplified the genomes using several SCoT molecular markers. We calculated of genetic diversity and polymorphism parameters using GenAlex 6.4, Geno-Dive ver.2, PopGene, PopART and Structure 2.3.4. The parameters of genetic polymorphism differed between the populations. Moreover, a low rate of gene flow supported a moderate level of population's genetic diversity and differentiation. According to haplotypes network (TCS) analysis, a high level of genetic mutation has occurred among the individuals of some populations leading to high intrapopulation diversity. On the basis of structure analysis and Nei's genetic distance, the examined populations were classified into four genotypes. CONCLUSIONS: The clustering pattern of the populations in each group was not related to geographical distance or phytogeography. It seems that the wide geographic distribution, a small gene flow rate and the occurrence of a high level of genetic mutation lead to infraspecific genetic differentiation in the species and we suppose some infraspecific ranks exist for it.

Genetic continuity of Indo-Iranian speakers since the Iron Age in southern Central Asia.

Since prehistoric times, southern Central Asia has been at the crossroads of the movement of people, culture, and goods. Today, the Central Asian populations are divided into two cultural and linguistic groups: the Indo-Iranian and the Turko-Mongolian groups. Previous genetic studies unveiled that migrations from East Asia contributed to the spread of Turko-Mongolian populations in Central Asia and the partial replacement of the Indo-Iranian populations. However, little is known about the origin of the latters. To shed light on this, we compare the genetic data on two current-day Indo-Iranian populations - Yaghnobis and Tajiks - with genome-wide data from published ancient individuals. The present Indo-Iranian populations from Central Asia display a strong genetic continuity with Iron Age samples from Turkmenistan and Tajikistan. We model Yaghnobis as a mixture of 93% Iron Age individual from Turkmenistan and 7% from Baikal. For the Tajiks, we observe a higher Baikal ancestry and an additional admixture event with a South Asian population. Our results, therefore, suggest that in addition to a complex history, Central Asia shows a remarkable genetic continuity since the Iron Age, with only limited gene flow.

Interspecific Gene Flow and Mitochondrial Genome Capture during the Radiation of Jamaican Anolis Lizards (Squamata; Iguanidae).

Gene flow and reticulation are increasingly recognized as important processes in the diversification of many taxonomic groups. With the increasing ease of collecting genomic data and the development of multispecies coalescent network approaches, such reticulations can be accounted for when inferring phylogeny and diversification. Caribbean Anolis lizards are a classic example of an adaptive radiation in which species have independently radiated on the islands of the Greater Antilles into the same ecomorph classes. Within the Jamaican radiation at least one species, Anolis opalinus, has been documented to be polyphyletic in its mitochondrial DNA, which could be the result of an ancient reticulation event or incomplete lineage sorting (ILS). Here, we generate mtDNA and genotyping-by-sequencing (GBS) data and implement gene tree, species tree, and multispecies coalescent network methods to infer the diversification of this group. Our mtDNA gene tree recovers the same relationships previously inferred for this group, which is strikingly different from the species tree inferred from our GBS data. Posterior predictive simulations suggest that our genomic data violate commonly adopted assumptions of the multispecies coalescent model (MSCM), so we use network approaches to infer phylogenetic relationships. The inferred network topology contains a reticulation event but does not explain the mtDNA polyphyly observed in this group; however, coalescent simulations suggest that the observed mtDNA topology is likely the result of past introgression. How common a signature of gene flow and reticulation is across the radiation of Anolis is unknown; however, the reticulation events that we demonstrate here may have allowed for adaptive evolution, as has been suggested in other, more recent, adaptive radiations. [Adaptive radiation; hybridization; introgression; multispecies network coalescent; posterior predictive simulation.].

Evaluation of gene flow and genetic diversity in rice accessions across Kurdistan region-iraq using SSR markers.

BACKGROUND: In recent years, farmers have complained that the only way to obtain seeds is to select plants that show good performance under local climate conditions in the region. This study aimed to investigate the diversity of rice accessions grown in the region to build a breeding program. METHODS AND RESULTS: A total of 62 accessions of rice from farmers and research stations were collected from the Kurdistan region, including short-grain and long-grain types, for molecular genetics and diversity analysis. In this study, 37 polymorphic simple sequence repeat (SSR) markers were selected with several molecular genetics software programs. The results show that these SSR markers are very effective for this investigation, generating a total of 152 observed alleles (Na), 75.166 Effective number of alleles (Ne) and an average of 4.1 and 2.03 alleles per locus, respectively. The average polymorphic information content (PIC) per locus was recorded as 0.404. The research presented here confirms two subpopulations, japonica (C1 and C2) and indica (C3), based on molecular genetics data analysis. Analysis of molecular variance revealed that the 72% variance was due to the variation among populations and 28% within the population. CONCLUSIONS: Altogether, these results indicate that there is very low gene flow. These results show the importance of the study of genetic diversity and relationships for starting breeding and improvement programs for rice in the Kurdistan region.

A genomic and morphometric analysis of alpine bumblebees: Ongoing reductions in tongue length but no clear genetic component.

Over the last six decades, populations of the bumblebees Bombus sylvicola and Bombus balteatus in Colorado have experienced decreases in tongue length, a trait important for plant-pollinator mutualisms. It has been hypothesized that this observation reflects selection resulting from shifts in floral composition under climate change. Here we used morphometrics and population genomics to determine whether morphological change is ongoing, investigate the genetic basis of morphological variation, and analyse population structure in these populations. We generated a genome assembly of B. balteatus. We then analysed whole-genome sequencing data and morphometric measurements of 580 samples of both species from seven high-altitude localities. Out of 281 samples originally identified as B. sylvicola, 67 formed a separate genetic cluster comprising a newly-discovered cryptic species ("incognitus"). However, an absence of genetic structure within species suggests that gene flow is common between mountains. We found a significant decrease in tongue length between bees collected between 2012-2014 and in 2017, indicating that morphological shifts are ongoing. We did not discover any genetic associations with tongue length, but a SNP related to production of a proteolytic digestive enzyme was implicated in body size variation. We identified evidence of covariance between kinship and both tongue length and body size, which is suggestive of a genetic component of these traits, although it is possible that shared environmental effects between colonies are responsible. Our results provide evidence for ongoing modification of a morphological trait important for pollination and indicate that this trait probably has a complex genetic and environmental basis.

Distinguishing gene flow between malaria parasite populations.

Measuring gene flow between malaria parasite populations in different geographic locations can provide strategic information for malaria control interventions. Multiple important questions pertaining to the design of such studies remain unanswered, limiting efforts to operationalize genomic surveillance tools for routine public health use. This report examines the use of population-level summaries of genetic divergence (FST) and relatedness (identity-by-descent) to distinguish levels of gene flow between malaria populations, focused on field-relevant questions about data size, sampling, and interpretability of observations from genomic surveillance studies. To do this, we use P. falciparum whole genome sequence data and simulated sequence data approximating malaria populations evolving under different current and historical epidemiological conditions. We employ mobile-phone associated mobility data to estimate parasite migration rates over different spatial scales and use this to inform our analysis. This analysis underscores the complementary nature of divergence- and relatedness-based metrics for distinguishing gene flow over different temporal and spatial scales and characterizes the data requirements for using these metrics in different contexts. Our results have implications for the design and implementation of malaria genomic surveillance studies.

Patterns of genetic diversity, spatial genetic structure and gene flow in Campomanesia xanthocarpa: insights from SSR markers of different genomic origins.

Campomanesia xanthocarpa (Mart.) O. Berg is a South American fruit tree species with important ecological and medicinal properties, which remnants are currently found mainly in isolated forest fragments. In this study, SSR markers from three different genomic origins (gene-linked, nuclear neutral, and organellar) were used to evaluate the patterns of genetic diversity, fine-scale spatial genetic structure and historical gene flow in fragmented forest formations of C. xanthocarpa from the Atlantic Forest in southern Brazil. Our results show that the forest fragments present moderate to high levels of genetic diversity in comparison to species presenting similar life traits, although a trend opposite to expected was observed concerning gene-linked and neutral SSR markers. The fine-scale spatial genetic structure revealed different patterns in short and large distance classes, with a distinct influence of gene-linked and neutral markers in driving the genetic structure in each distance class. The presence of an isolation-by-adaptation pattern implies the need for maintenance of the current remnants to assure the conservation of the private alleles. Finally, as the genetic diversity is found predominantly within forest fragments, programs of seed collection and/or genetic rescue should prioritize a larger number of individuals within each fragment, to increase the sampled diversity.

Against the Odds: Hybrid Zones between Mangrove Killifish Species with Different Mating Systems.

Different mating systems are expected to affect the extent and direction of hybridization. Due to the different levels of sexual conflict, the weak inbreeder/strong outbreeder (WISO) hypothesis predicts that gametes from self-incompatible (SI) species should outcompete gametes from self-compatible (SC) ones. However, other factors such as timing of selfing and unilateral incompatibilities may also play a role on the direction of hybridization. In addition, differential mating opportunities provided by different mating systems are also expected to affect the direction of introgression in hybrid zones involving outcrossers and selfers. Here, we explored these hypotheses with a unique case of recent hybridization between two mangrove killifish species with different mating systems, Kryptolebias ocellatus (obligately outcrossing) and K. hermaphroditus (predominantly self-fertilizing) in two hybrid zones in southeast Brazil. Hybridization rates were relatively high (~20%), representing the first example of natural hybridization between species with different mating systems in vertebrates. All F1 individuals were sired by the selfing species. Backcrossing was small, but mostly asymmetrical with the SI parental species, suggesting pattern commonly observed in plant hybrid zones with different mating systems. Our findings shed light on how contrasting mating systems may affect the direction and extent of gene flow between sympatric species, ultimately affecting the evolution and maintenance of hybrid zones.

Deforestation is the turning point for the spreading of a weedy epiphyte: an IBM approach.

The rapid spread of many weeds into intensely disturbed landscapes is boosted by clonal growth and self-fertilization strategies, which conversely increases the genetic structure of populations. Here, we use empirical and modeling approaches to evaluate the spreading dynamics of Tillandsia recurvata (L.) L. populations, a common epiphytic weed with self-reproduction and clonal growth widespread in dry forests and deforested landscapes in the American continent. We introduce the TRec model, an individual-based approach to simulate the spreading of T. recurvata over time and across landscapes subjected to abrupt changes in tree density with the parameters adjusted according to the empirical genetic data based on microsatellites genotypes. Simulations with this model showed that the strong spatial genetic structure observed from empirical data in T. recurvata can be explained by a rapid increase in abundance and gene flow followed by stabilization after ca. 25 years. TRec model's results also indicate that deforestation is a turning point for the rapid increase in both individual abundance and gene flow among T. recurvata subpopulations occurring in formerly dense forests. Active reforestation can, in turn, reverse such a scenario, although with a milder intensity. The genetic-based study suggests that anthropogenic changes in landscapes may strongly affect the population dynamics of species with 'weedy' traits.

And yet, migration, population growth or mortality are not balanced among Cantabrian brown bear subpopulations. Reply to Blanco et al (2020).

In a recent paper, we presented new evidence and provided new insights on the status of Cantabrian brown bear subpopulations, relevant for this species conservation. Namely, we revealed the likely phylogeographic relation between eastern Cantabrian subpopulation and the historical Pyrenean population. We have also detected an asymmetric flow of alleles and individuals from the eastern to the western subpopulation, including seven first-generation male migrants. Based on our results and on those of previous studies, we called the attention to the fact that Eastern Cantabrian brown bears might be taking advantage of increased connectivity to avoid higher human pressure and direct persecution in the areas occupied by the eastern Cantabrian subpopulation. In reply, Blanco et al (2020) [11] have criticized our ecological interpretation of the data presented in our paper. Namely, Blanco and co-authors criticize: (1) the use of the exodus concept in the title and discussion of the paper; (2) the apparent contradiction with source-sink theory; (3) the apparent overlooking of historical demographic data on Cantabrian brown bear and the use of the expression of population decline when referring to eastern subpopulation. Rather than contradicting the long and growing body of knowledge on the two brown bear subpopulations, the results presented in our paper allow a new perspective on the causes of the distinct pace of population growth of the two brown bear subpopulations in the last decades. Here, we reply to the criticisms by: clarifying our ecological interpretation of the results; refocusing the discussion on how the new genetic data suggest that currently, the flow of individuals and alleles is stronger westward, and how it may be linked to direct persecution and killing of brown bears. We provide detailed data on brown bear mortality in the Cantabrian Mountains and show that neither migration, gene flow, population increase nor mortality are balanced among the two subpopulations.